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Abstract 



Quantum theory of interference phenomena does not take the di- 
ameter of the particle into account, since particles were much smaller 
than the width of the slits in early observations. In recent experiments 
with large molecules, the diameter of the particle has approached the 
width of the slits. Therefore, analytical description of these cases 
should include a finite particle size. The generic quantum interference 
setup is an asymmetric double slit interferometer. We evaluate the 
wave function of the particle transverse motion using two forms of the 
solution of Schrodinger's equation in an asymmetric interferometer: 
the Fresnel-Kirchhoff form and the form derived from the transverse 
wave function in the momentum representation. The transverse mo- 
mentum distribution is independent of the distance from the slits, 
while the space distribution strongly depends on this distance. Based 
on the transverse momentum distribution we determined the space 
distribution of particles behind the slits. We will present two cases: 
a) when the diameter of the particle may be neglected with respect 
to the width of both slits, and b) when the diameter of the particle is 
larger than the width of the smaller slit. 
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1 Introduction 



Until recently various quantum interference experiments were conducted 
with objects (photons, electrons, neutrons,...) of the size much smaller than 
the characteristic dimensions of the diffraction structure [1]. The first single 
slit experiment with Rydberg atoms, objects of non negligible size with re- 
spect to the width ot the slit was performed by Fabre et al. [2]. By measuring 
transmission through micrometer size slits, these authors determined the size 
of Rydberg atoms. Later, Hunter and Wadlinger [3] proposed the single-slit 
diffraction experiment in order to measure the diameter of the photon. 

In order to investigate the reasons of unobservability of quantum effects 
in the classical world Arndt et al. [4], Nairz et al. [5], Brezger et al. [6] 
performed quantum interference experiments with objects of large mass and 
diameter, including macromolecules. The experiments raise various questions 
about theoretical concepts. The following three are obvious: 

1) Are there new effects applying to particles bigger than their de Broglie 
wavelength? 

2) Does internal structure have influence on interference? 

Arndt et. al. emphasized the task [4]: "Here we report the observation 
of de Broglie wave interference of Ceo molecules by diffraction at a material 
absorption grating. This molecule is the most massive and complex object in 
which wave behavior has been observed. Of particular interest is the fact that 
Ceo is almost a classical body, because of its many excited internal degrees 
of freedom...". 

3) What happened with an ensemble of incoming particles if slits are 
smaller than the diameter of the particles? 

In the experiment of Arndt et. al. the de Broglie wavelength of the inter- 
fering fullerenes is already smaller than their diameter by a factor of almost 
400 and authors pointed out that "it would be certainly interesting to in- 
vestigate the interference of objects the size of which is equal or even bigger 
than the diffracting structure" [4]. 

These experiments could shed more light on a long standing dilema whe- 
ther each quanton consists of a particle and accompanied wave (as two differ- 
ent compatible entities) [7], or quantons sometimes behave like a wave and 
sometimes behave like a particle (obeying principle of complementarity) [8]? 
The following citations illustrate the present situation. 
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"We performed an experiment which was proposed by Ghose, Home and 
Agarwal showing both classical wave-like and particle-like behaviors of single 
photon states of light in a single experiment, in conformity with quantum 
optics." [9]. "Here we give a detailed justification of our claim that this 
experimental results contradict the tenet of mutual exclusiveness of classical 
wave and particle pictures assumed in Bohr's complementarity principle." 
[10]. "Simultaneous observations of wave and particle behavior is prohibited" 

:n]- 

"Although interference patterns were once thought of as evidence for wave 
motion, when looked at in detail it can be seen that the electron arrive in 
individual lumps. ... We must therefore conclude that electrons show wave- 
like interference in their arrival pattern despite the fact that they arrive in 
lumps just like bullets". [12] 

"It is frequently said or implied that the wave-particle duality of matter 
embodies the notion that a particle - the electron, for example - propagates 
like a wave, but registers at a detector like a particle. Here one must again 
exercise care in expression so that what is already intrinsically difficult to 
understand is not made more so by semantic confusion. The manifestations 
of wave-like behavior are statistical in nature and always emerge from the 
collective outcome of many electron events... That electrons behave singly 
as particles and collectively as waves is indeed mysterious, ... " [13] 

"Each atom is therefore at the same time a particle and a wave, the wave 
allowing one to get the probability to observe the particle at a given place." 
[14] 

"...Ever since then the two sides of the same quantum object appeared 
together: on the one hand the non-local wave nature needed to describe the 
unperturbed propagation and on the other hand the local aspect of the object 
when it is registered by the detector" [4] 

In this paper we present the theoretical study of the dependence of the 
quantum interference pattern on the diameter of the particle, assuming that 
the characteristic sizes of the diffraction structure are of the order of the 
diameter of the particle. We argue that an asymmetric double-slit interfer- 
ometer (an interferometer whose slits have different widths 5i and 62 ) is the 
generic case for this study. 
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2 The particle wave function behind an 
asymmetric grating 



We shall now determine, the wave function of a quanton which travels 

— * 

with velocity, v — vi through the region I, towards the slits and is then 
sent through the slits to the region II (Fig. 1). Results in this section are 
valid for arbitrary slits. This wave function is a stationary solution of the 
time-dependent two dimensional Schrodinger equation 

h 2 ( d 2 d 2 \ d 

(a? + ^J* feB - f) = ^ (1) 

The solution of (1) has the form 

*(x,y,t) = e- iut <p(x,y), (2) 

where Tiuj = mv 2 /2 and p = mv = Kk. Space dependent function ip(x,y) 
satisfies the Helmholtz equation 

h 2 ( d 2 d 2 \ 
'^{d^ + dy- 2 ) ViX,y) = hU ^ y) - (3) 

The solution of this equation in the region I is a spherical wave 

V (P')= V (x\y')=A e ^ ri (4) 

where A is a constant and r' is a distance (Fig. 1) from the source (P ) to 
the point P' = (x', y') in the region I. The distance a of the double-slit screen 
from the source Pq being very large compared to the width of the slits, this 
spherical wave at the slit points (x' = x",y' = 0) may be approximated by 
the plane wave. In the region II the equation (3) is as simple as before but 
initial condition makes the solution more difficult. 

Solution known as Fresnel-Kirchhoff diffraction formula [15] reads: 



V(x,y) = -— dx" [1 + cosx], (5) 

IX a J a s 



ika „ „iks 
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where s = yy 2 + (x" — x) 2 , cosx = y/s, A = 2ix/k. The region A is the 
union of all intervals along the x-axis where slits are open. From now on x" 
represents a variable of integration along the line of the slits. 

Far enough from the slits wave function resembles the Fourier transform 
of the wave field accross the aperture. This can be verified from the Fresnel- 
Kirchhoff solution. In the far region it follows: 

A p ika iky ,. 

<p( x ,y) » J± e — e — / dx"(p(x",0)e- ikxx l y . (6) 
A a y J a 

Wave function is now separable into two functions, one depending on y and 
the other depending on K x = kx/y [16]: 

<p(x,y)=D(y)T(K x ) 

iAe ika e iky 

D{y) = -v^F-r 

A a y 

t{k x ) = 4rl dx " ^ x "> °) e ' tKxX "- ( ? ) 

V 2ir J A 

The solution of equation (3) in region II can be written in another form 
[17, 18, 19]. This form is more convenient for our analysis than the form (5). 
With approximation valid for small diffraction angles x we have: 

tp(x, y) = e iky -= / dk x c(k x ) e^e-^ = e^^x, y). (8) 

V Z7T J -co 



. 2tt 

where c(k x ) is the Fourier transform of the function tp(x,y) on the aperture 
<p(x,y = 0): 

j r+co 

c{k x ) = -= / dx" <p(x", 0) e~ ik * x . (9) 



2,TT 

Inserting (9) into (8), after integration over k x one finds 

<p(x, y) = e ~^ etky ^^j j A Viz", ^S)e kJ ^^dx". (10) 

This function is normalized \(p(x" ,y)\ 2 dx" = 1, provided J A \(p(x", 0)| 2 
dx" = 1. The form (10) clearly expresses wave function's dependence on the 
boundary condition and it is appropriate for numerical computation. 
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For large values of y the function ip(x,y) in (10) is approximated by 



<p( x ,y) = M e-ie^V^ 2 /^ / <p(x",0)e- ikxx "/ y dx". (11) 
V 2iry J a 

Taking Eq. (9) into account, Eq. (11) takes the form 

<p(x, y) = e iky ^e-^e ikx2/2y c(kx/y). (12) 

We see that the variable K x — y = ^ plays the role of k x . 

Since K x is proportional to x/y functions \<p(x, y) = const\ are family of 
functions of x spreading along the x— axis as y increases. In fact, for each 
value of \tp\, in the far field there exists the straight line with origin at the 
center of the grating along which this particular value of \ip\ propagates. 

3 The understanding of the space 

distribution using transverse momentum 
distribution 



By assuming that the motion of an atom along the y-axis may be 
treated classically and that the transverse motion is quantum, one is tempted 
to use the relation y = vt and to determine the time dependent function of 
the transverse motion ip(x,t) from the function <f)(x,y), by the following 
definition: 

if>(x, t) = <j>{x, vt) = f dk x c(k x )e lkxX e~ luJxt (13) 

where u x = hk 2 x /2m. We see that the function ip(x, t) has the form of a gen- 
eral solution of the one-dimensional Schrodinger equation. The wave function 
c(k x ) is then seen as a wave function of this one- dimensional (transverse) mo- 
tion in the momentum representation. It's modulus square, |c(/ca;)| 2 , deter- 
mines the distribution of transverse momenta. The wave function ^(x,y,t) 
from Eq. (2) is expressed through vp(x,t) as 

V(x, y, t) = e iky e~ iu}t ^(x, t). (14) 
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By taking Eq. (12) into account one concludes that in the far field the 
relation (13) between the wave functions ip(x,t) and c(k x ) reduces to the 
simplier form: 

V> (*, t = ^) = ^-e-^e^ 2 ^c{kx/y). (15) 

Based on the above factorization of the wave function ^(x, y, t) and the prop- 
erties of its factors summarized above, we proposed [18] the new expression 
for the probability density P(x,t) fo the particle's arrival to a certain point 
(x, y = vt) at time t: 

P (x,^j = P(x,t) EE 

= f °° dk x f °° dx" \c n {k x )\ 2 \<P{x\ 0)| 2 <5 (x - x" - ^] . (16) 

Particles emerge from different points (x", 0) at the aperture. That is the 
reason for integration over x". The contribution of each point at the aperture 
is proportional to |0(x",O)| 2 . The integration over dk x and the function 
l c n(&z)| 2 reflect the contribution of various angles/momenta in diffraction. 
Finally, 5-function assumes straight trajectory from a point (x" , 0) at the 
slits to the point (x, y) and leads to the simplified form 

r+oo 

P(x,t) = / dk x \c n (k x )\ 2 

J — CO 

By assuming that the function <f>(x", 0) = for x" £ A and <f>(x", 0) = 
const such that J A \4>(x", 0)| 2 = 1, for x" e A, the Eq. (17) is transformed to 
the following usefull form 

P(x,t)= Yl I 1 ^|c(**)| 2 = ^Pi(x,t). (18) 

Here x\ and the coordinates of the left and right edge of the i-th slit. 

The total probability density P(x, t) is a sum of n terms, Pi(x, t). Pi(x, t) 
is interpreted to be the probability that a quanton reaches (x, y = vt) at time 
t after passing through the i-th slit of the n-slits grating. 



X 



m 



-,o 



(17) 
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Numerical calculation shows that far from the slits the function P(x, t) 
(Fig. 4) is very nearly equal to the exact probability density \^(x, y, t)\ 2 = 
\ip(x, t)\ 2 (Fig. 2). Near the slits P(x,t) and \ip(x,t)\ 2 qualitatively look 
similarly but they differ numerically. 



4 On the possible influence of particle's 
diameter on the interference pattern 



We outline an approach to investigate how the widths of the slits 
influence the interference pattern in the double-slit experiment with quantons 
- photons, electrons, neutrons, atoms, molecules. 

Interference effects are visible when the wavelength of quantons is of 
the order of the distance between the slits. In practice, this distance is 
d — (2 — 50) A. The slit width is often equal or up to ten times smaller 
than the distance between the slits. In quantum interference experiments 
with electrons and neutrons the diameter of the particle is smaller than the 
wavelength. Consequently, in classical experiments the width 5 of the slits 
is much greater than D. But, depending on the velocity, atoms may have 
de Broglie wavelength which is smaller than the diameter of the atom. With 
macromolecules such a situation encounters more often, as shown in the ex- 
periment of Arndt et. al. and discussed by Arndt et al. [4] and Nairz et al. 
[5]. So, interference experiments with such quantum particles could have the 
slit widths smaller than the particle diameter. 

This requires a theoretical approach to quantum interference which takes 
the diameter of the particle into account [20]. A study of a quantum particle 
in an asymmetric double-slit interferometer (5i > 5 2 ) seems to be useful for 
this purpose because we identify two characteristic cases for the ratio of slit 
widths 5i and 5 2 and the diameter of the particle D: 

a) The diameter D is negligeable with respect to the widths 5\ and 5 2 ■ 

b) The diameter D is greater than the width 5 2 , where 5 2 < 5±. 

In the case a), which was until recently the only case of physical interest, 
there is no need to consider or take into account the diameter of the particle. 
The particle momentum |c(/c x )| 2 and space distribution \ip(x, t)\ 2 behind the 
grating are determined by the wave function 




(19) 
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where 



i 



0(x,O) = 4>(x,0) 







x e A, 
x g A 



■d-Si 



d-S 2 



and 



c(k x 



^/2tt(5i + 6 2 ) k * . 



e lk * d/2 sin 



Mi _|_ e ~ik x d/2 gin ^2 



(20) 



(21) 



The functions \ip(x, t)\ 2 and |c(/c x )| 2 are graphically represented at Fig. 2 and 
Fig. 3, for the chosen set of parameters. 

In the case b), we are faced with the question how and where to take 
the diameter of the particle into account. We know that the diameter of the 
particle is not incorporated anywhere in the Schrodinger equation. But, we 
expect that a particle with diameter D, such that 5i > D > 5 2 could not 
pass through the second slit. 

So, it seems to us that we are forced to assume that wave functions in 
the coordinate and momentum representation in the case b) is also given by 
expressions (19)-(21). 

The momentum distribution |c(/c x )| 2 of particles is given also by (21), 
because it is determined by the values of the wave function at the boundary. 

But the space distribution of particles in case b) is different from the 
space distribution in case a), because the particles arriving to the smaller 
slit can not go through. We conclude that particle distribution in case b) is 
given by Pi(x,t) from the expression (18) of P(x,t). 



P(x,t) w AOm) 



1 



f (*-*}) 



(22) 



The probability Pi(x,t) is graphically represented in Fig. 5. 

5 Conclusion 



Inspired by current efforts to perform diffraction and interference experiments 
with objects of size that is equal or even larger than the diffraction structure, 
we outline an approach to investigate how the particle diameter influences 
the interference pattern in an asymmetric double slit interferometer. 
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We identify two characteristic cases for the ratio of slit widths 5\ and 82 
and the diameter D of the particle: a) D <C Si and D <C 62, b) 5\ > D > 
5 2 - The wave function behind the grating has the same form in both cases 
because it is the solution of the Schrodinger equation which is not sensitive 
to the diameter of the particle. 

The space distribution of particles in case a) is given as usual by the 
modulus square of this function. Using the same wave function and assuming 
that a particle with diameter D, such that 5\ > D > 62 could not pass 
through the second slit, we determine the space distribution in case b). We 
conclude that the momentum distribution of particles behind the grating is 
the same in cases a) and b). 
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Figure captions 



Fig. 1. Illustration of a grating with n slits of various widths. 

Fig. 2. The particle distribution function \ip(x, t)\ 2 behind the asymmetric 
double slit grating (Si = 1 /im, S2 = 0.25 /im, d = 8 /im) close to the slits 
(a, b) and far from the slits (c, d). It is evaluated from the form (19) of the 
wave function. The initial longitudinal wave vector is k = 4n ■ 10 10 m -1 , the 
particle mass is m = 3.8189 ■ 10~ 26 kg. 

Fig. 3. The particle transverse momentum distribution |c(A; x )| 2 behind 
the asymmetric double-slit grating (Si = 1 /im, S2 = 0.25 /im, d = 8 /im). 

Fig. 4. The probability density P(x, t) of particles arrival to the point 
x at time t (y = vt) behind the asymmetric double slit grating (<5i = 1 /im, 
5 2 = 0.25 /im, d = 8 /im) close to the slits (a, b) and far from the slits (c, d). 
It is evaluated from Eq. (18). Particles' diameter D is negligible with respect 
to the widths of the slits. The initial longitudinal wave vector is k = 4tt ■ 
10 10 m _1 , the particle mass is m = 3.8189 • 10~ 26 kg. 

Fig. 5. The probability density Pi(x, t) of particles reaching (x, y) at time 
t after passing through the larger slit, near the slits (a, b) and far from the 
slits (c, d). It is evaluated from Eq. (22). D is assumed to be larger than 5 2 
and smaller than 5±. The values of parameters are the same as in captions 
of Figs. 2,3,4. 
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